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We calculated the energy spectra and the fluxes of electron neutrinos in extragalactic space 
emitted in the process of the evaporation of primordial black holes (PBHs) in the early universe. It 
was assumed that PBHs are formed by a blue power- law spectrum of primordial density fluctuations. 
In the calculations of neutrino spectra the spectral index of density fluctuations and the reheating 
temperature were used as free parameters. The absorption of neutrinos during propagation in the 
space was taken into account. We obtained the bounds on the spectral index assuming validity of 
the standard picture of gravitational collapse and using the available data of several experiments 
with atmospheric and solar neutrinos. The comparison of our results with the previous constraints 
(which had been obtained using diffuse photon background data) shows that such bounds are quite 
' sensitive to an assumed form of the initial PBH mass function. 

o 

I. INTRODUCTION 

b 

, Some recent inflation models (e.g., the hybrid inflationary scenario predict the "blue" power-spectrum of 
' primordial density fluctuations. In turn, as is weU known, the significant abundance of primordial black holes (PBHs) 
is possible just in the case when the density fluctuations have an n > 1 spectrum (n is the spectral index of the initial 
''si" ■ density fluctuations, n> 1 spectrum is, by definition, the "blue perturbation spectrum"). 

^ ', Particle emission from PBHs due to the evaporation process predicted by Hawking |^ may lead to observable 
effects. Up to now, PBHs have not been detected, so the observations have set limits on the initial PBH abundance 
or on characteristics of a spectrum of the primordial density fluctuations. In particular, PBH evaporations contribute 
r to the extragalactic neutrino background. The constraints on an intensity of this background (and, correspondingly, 

■ on an PBH abundance) can be obtained from the existing experiments with atmospheric and solar neutrinos. The 
obtaining of such constraints is a main task of the present paper. 

■ The spectrum and the intensity of the evaporated neutrinos depend heavily on the PBH's mass. Therefore, the 
\ great attention should be paid to the calculation of the initial mass spectrum of PBHs. We use in this paper the 

rS i' following assumptions leading to a prediction of the PBH's mass spectrum. 

1. The formation of PBHs begins only after an inflation phase when the universe returns to the ordinary radiation- 
dominated era. The reheating process is such that an equation of state of the universe changes almost instantaneously 
^ ' into the radiation type (e.g., due to the parametric resonance JH] ) after the inflation. 

, 2. It is assumed, in accordance with analytic calculations [Epj that a critical size of the density contrast needed for 
the PBH formation, 5c, is about 1/3. Further, it is assumed that all PBHs have mass roughly equal to the horizon 
mass at a moment of the formation, independently of the perturbation size. 

3. Summation over all epochs of the PBH formation can be done using the Press-Schechter formalism tel. This 



O 



X 



5h \ formalism is widely used in the standard hierarchial model of the structure formation for calculations of the mass 
- - I distribution functions (see, e.g., 0). 

It was shown recently that near the threshold of a black hole formation the gravitational collapse behaves as a 
critical phenomenon . In this case the initial mass function will be quite different from that which follows from the 
standard calculations of refs. |^,0.The first calculations of the PBH initial mass function (e.g., in ref. [^) have been 
done under the assumption that all PBHs form at the same horizon mass (by other words, that all PBHs form at 
the smallest horizon scale immediately after reheating). The initial PBH mass spectrum for the case of the critical 
collapse, based on the Press-Schechter formalism, was obtained in the recent work 

The calculations in the present paper are based on the standard [^|j5j picture of the gravitational collapse leading 
to a PBH formation. The case of the critical collapse will be considered in a separate work. 
The plan of the paper is as follows. 

In Sec.^ we give, for completeness, the brief derivation of a general formula for the initial PBH mass spectrum. 
The final expression is presented in a form, which is valid for an arbitrary relation between three physical values: 
the initial PBH mass, the fluctuation mass (i.e., the mass of the perturbed region) at a moment of the collapse, 
and the density contrast in the perturbed region (also at a moment of the collapse). This expression contains the 
corresponding results obtained in refs. [^|J|] as particular cases. As in refs. 0,^, the derivation is based on the linear 
perturbation theory and on the assumption that a power spectrum of the primordial fluctuations can be described by 
a power law. 
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In Sec, II] we derive the approximate formula for a calculation of the extragalactic neutrino background from PBH 
evaporations. We do not use cosmological models of the inflation and of the spectrum of primordial fluctuations, so 
there are two free parameters: the reheating temperature and the spectral index. At the end of the section, some 
examples of instantaneous neutrino spectra from the evaporation of an individual black hole are presented. 

In Sec, IV we perform the normalization of the primordial power spectrum of density perturbations by COBE data 
on the angular power spectrum of CMB temperature fluctuations. We consider the case of flat cosmological models 
with nonzero cosmological constant. We show, that the result of this normalization rather weakly depends on the 
matter content of the present Universe, Q,m = 1 — ^a- 

In Sec.^ we present the results of numerical calculations of the neutrino background spectra from PBH evaporations, 
with taking into account the effects of neutrino absorption in the space. For background neutrino energies '--^ 10 — 
100 M eV we studied relative contributions to the background intensity of different cosmological rcdshifts (and it is 
shown that at high reheating temperatures the characteristic values of the redshifts are very large) . 

In Sec. VI we consider the possibilities of constraining the spectral index, using experiments with neutrinos of natural 
origin. 



In Sec. VII the spectral index constraints, followed from our calculations and from available data of the neutrino 
experiments, are given. 



II. THE INITIAL MASS SPECTRUM OF PBHS 

As is pointed out in the Introduction, we use the Press-Schechter formalism which allows to carry out the 
summation over all epochs of PBHs formation. According to this formalism, the mass spectrum of density fluctuations 
( i.e., the number density of regions with mass between M and M + dM ) is calculated by the formula 



n (M, 5^) dM 
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dM 



Here, P{M, Sc) is the fraction of regions having sizes larger than R and density contrast larger than Sc 
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In these formulas, 5 is the initial density contrast, ue. is the stan dar d deviation of the density contrast of the regions 
with a comoving size R and mass M. The factor of 2 in Eq.( |2.2| ) ensures the correct normalization of n(M, i5c), 
integrated over all M . This factor is usually added to take into account (approximately) the contribution of underdense 
regions to gravitationally bound objects 0|. 

It is convenient to introduce the double differential distribution n(M, 5), the integral over which gives the total 
number of fluctuated regions: 



n{M, 6)dMdS 

Using Eqs.( ^.l| - |2.3[ ), for n{M,6) one has the expression 

' 2 1 dan 
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To obtain the mass spectrum of PBHs one must introduce the Mbh variable. Besides, we will use the variable iJ', 
connected with 5 by the relation 



2/3 



(2.6) 
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Here, Mi is the horizon mass at the moment of a beginning of the growth of density fluctuations. This new variable 
is the density contrast at the moment of the collapse. The new distribution function is 



nBH(MBH,0 ) = n(MBH,0) — - 



dS' dMsH 

Further, we assume that there is some functional connection between Mbh, M and S': 

MBH^f{M,d') . 



In this case one can rewrite Eq.(2.7) in the form 
tibh^Mbh, 5') = ni M, f^' 



-2/3N 



M 
M 



-2/3 



df{M, S')/dM 



and the PBH mass spectrum is given by the integral 

riBHiMBH) = j nBH{MBH,S')dS' . 

Now, to connect the mass spectrum with the spectral index one can use the relations 

MX -2/3 



(yR = cfH{M) 
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In these relations aniM) is the horizon crossing amplitude, i.e., the standard deviation of the density contrast in 
the perturbed region at the moment when this fluctuation crosses horizon. The derivation of Eqs.(2.11) is given, for 
completeness , in Sec JTVI. 



Using Eqs.(2.11) in the expression (2.5) for n{M,5) one obtains 



^ + 3 /2 ft 1 f M 



-4/3 



Substituting this expression in the r.h.s. of the Eq.(p.9|) and using Eq( |2.10|) one has 



nBH[AlBH) = „ • \ -Pi 
6 V TT 



M2 



'H 
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Values of M in r.h.s of Eq.(2.13) are expressed through Mbh, 5' by the relation (| 

Formula ( 2.13 ) is the final expression for the PBH mass spectrum, the main result of this Section. It is valid for 
any relation between PBH mass Mbh, mass of the original overdense region M and the density contrast 5' . 

In the particular case of a near critical collapse 



therefore 



/(M, 5') = Ml'^M'^'^k{5' - 5^y^ = CM^/^M2/3 



M = A4/^A/ri/2^-3/2 ; df /dM ^ -Ml^'^Mgll'^i^/'^ 



(2.14) 



(2.15) 



Eq.(2.14) can be rewritten in the form, derived in the ref. using the connection between M and the horizon 
mass Mh (which is equal to a fluctuation mass at the moment when the fluctuation crosses horizon): 



Mh = M^^M^/^ 



The resulting formula for the spectrum 

n + 3 12 
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coincides, as can be easily proved, with the expression derived in ref. 

The mass spectrum formula for the Carr-Hawking collapse can be obtained analogously, with the relation /(Af, S') 
^i/2jvf V3^2/3 (ggg Eq.( ^.12|) below), or from Eq. (|2.17|) using the substitutions 



7fc 



, k 



vl/2 



and the approximate relation 



dS' 
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In this case one obtains the expression derived in ref. 

nBH{MBH) = 



+ 3 /2 ,,1/2,, -5/2 -1 
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(2.18) 
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One can see from Eqs.(2.17) and ( 2.201 ) that the PBH mass spectrum following from the Press-Shechter formalism 
has quasi power form in both considered cases ^^bh'^^ ■ ^'^ Carr-Hawking case M'j!^}p ~ Mh , in contrast with 



this in the critical collapse case it is possible that Mbh ^ Mh 



III. FORMULA FOR NEUTRINO DIFFUSE BACKGROUND FROM PBHS 



The starting formula for a calculation of the cosmological background from PBH evaporations is ||1 
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S{E)= I n,o^^^—^f{E{l + z))dV,om ■ (3.1) 



Here , ricom is the comoving number density of the sources (in our case the source is an evaporating PBH of the 
definite mass m ), oq is the scale factor at present time, t = to, f{E) is a differential energy spectrum of the source 
radiation. Vcom is a comoving volume of the space filled by sources, therefore 

dVcom. = Oq / =dfl . (3.2) 

Here, k is the curvature coefficient, and p is the radial comoving coordinate. Using the change of the variable, 

dp dt 



the comoving number density can be expressed via the initial density Ui 



(3.3) 



ito)=nJ^y(^Y = nJ'-^y . (3.4) 
\ a / \an J Van/ 



^com — '^phys 



ao J V ao 



Substituting Eqs. (3.2) - (3.4) in Eq. (3.1) one obtains 

3 

a \ flo 

In our concrete case the source of the radiation is a Hawking evaporation 



S{E)^n, I dt^i^] fiE{l + z)) . (3.5) 



n,f{E{l + z))^JdmnBHim,t)fHiE{l + z),m) . (3.6) 
Here , UBHi'tn, t) is the PBH mass spectrum at any moment of time, fniE, m) is the Hawking function 

^"^'''"^^^^^^^^(l^)-(-l)- ■ ^'-'^ 
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Here, E is the energy of an evaporated particle (it lies in the interval [0,oo] for massless particles), rs{E,m) is the 
coefBcient of the absorption by a black hole of a mass m, for an particle having spin s and energy E. 

In our calculations (we considered only the simplest case of a nonrotating uncharged hole) we used the approximation 
when Ts depends only on mE (as is the case for massless particles), i.e. we ignore the decrease of Tg at nonrelativistic 
energies of evaporated particles. But, working in this massless limit, we nevertheless used in a low mE region the exact 
forms of Fg, which are obtained numerically (concretely, we used the plots for absorption cross section, a « Ts/E"^, 
as a function of mE, for s ~ 1/2 and s = 1, from the Fig.l of work |l^), rather than the low-energy limits obtained 
analytically in works 

The initial spectrum of PBHs is given by Eq. (2.20). The minimum value of PBH mass in the specrtum can be 
obtained from the following considerations. It is known that the condition of a recoUapse of the perturbed region 
is the following: radius of a maximum expansion of this region, Rc, must obey the inequality 



R, 



(3.8) 



Here, Rh is the horizon radius at ic, the moment of time, when the radius of the perturbation region is equal to 
maximum one. It is known also [p[ that tc and the initial moment of time, are connected by the relation 



(3.9) 



where 5i is the initial density contrast for which there is the condition following from Eq.(3 



-2/3 



(3.10) 



Using Eq.(3.8) and the formula for horizon mass at 



one obtains the connection between Mbh and t^: 



Mbh = ^^Rlp{U)^-f^'HU^-^ 



The use of Eqs.( ^!9plo| ) in Eq.( |3.12| ) leads to the relation 

Mbh = ^^'^mI'Hi^/^ 



(3.11) 



(3.12) 



(3.13) 



From Eq.(^.13D one follows that the minimum value of a PBH mass in the PBH mass spectrum is given by the simple 
formula 



(3.14) 



Strictly speaking, the minimum value of |^ is equal to 1 (corresponding to (5™"^ — 1). Therefore, the absolute 
minimum value of PBH mass, following from Eq.(3.12), is equal to However, the use of such a value as a 



border value of the PBH mass spectrum is clearly inconsistent. It contradicts with the relation ( p. IS ) and, in general, 
with the fact that, due to the ex ponen tial damping Q, the bulk of the PBH intensity is determined just by the 
minimum value of 6i (given by Eq.( 3.10| )). 

To take into account the existence of the minimum we must add to the initial spectrum expression the step factor 
^{Mbh — "f^^'^Mi). The connection of the initial mass value Mbh and the value at any moment t is determined by 
the solution of the equation p2[ 



dm 



a{m) 



(3.15) 



The function a(m) accounts for the degrees of freedom of evap orated particles and determines the lifetime of a black 
hole. In the approximation a = const the solution of Eq.( |3.15| ) is 



Mbh = (3at 



,1/3 



(3.16) 
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This decrease of PBH mass leads to the corresponding evohition of a form of the PBH mass spectrum. At any moment 
one has 



?^s_ff (™, t)dm 



{Sat + m^) 



3^2/3 



Sat 



1/3 



dm. 



(3.17) 



Substituting Eqs. (3.7), (3.17) in the integral in Eq.(3.6) , we obtain the final expression for the spectrum of the 
background radiation: 



S{E) = 



47r 



dt 



ao I ai 
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(3.18) 



(m - ((7'/'M,)3 - 3at) 



gufH{E{l + z),m) 



Here, g^, is a number of degrees of freedom. In the following we will interest in {ve + i^e) - spectrum, so gi, = 2. 
One should note that the corresponding expressions for the spectrum in refs. |]l^,^ contain the factor instead 

of the correct factor (^f^^ ■ It leads to a strong overestimation of large z contributions in S{E) (see below, Fig.|^). 
It is convenient to use in Eq.( |3.18 ) the variable z instead of t. In the case of flat models with nonzero cosmological 



constant one has 



dt 
dz 



Ho{l + z) 



(r!™(l + z)3 + a.(l + 2)^ + f^A) 



-1/2 



(3.19) 



Vlr = (2.4 • lO^h^Y^ , Q.„, = 1 - Vl,. - flA 



(3.20) 



The factor (f^)"^ can be expressed through the value of t^q, the moment of matter-radiation density equality: 



(1 + Zeq) 



-3 



3/2 



-3/2 



4f,2\-3/4 3/2 



The constant Hq, entering the Eqs.( |3.19y3.2l| ) is connected with h by definition: 

Ho = lOO/i km s"^ Mpc"^ . 



(3.21) 



(3.22) 



Throughout the paper we suppose that h = 0.67. One can see from Eqs. ( 3. 21 , 3.22 ) that the factor (^f^j does not 

depend on h and flm, as it must be. 

Integrating over PBH's mass in Eq.( 3.18| ), one obtains finally, after the change of the variable < on z, the integral 
over z: 



S{E) - J F{E,z)dlog^oiz + l) 



(3.23) 



In analogous calculations of the photon diffuse background integral over z in the expression for S{E) is cut off at 
z = zo ~ 700 because for larger z the photon optical depth will be larger than unity |l^ ]. In contrast with this , 
interactions of neutrinos with the matter can be neglected up to very high values of z. Therefore the neutrino diffuse 
background from PBH evaporations is much more abundant. The neutrino absorption effects are estimated below, in 
Sec.0. 

The evaporation process of a black hole with not too small initial mass is almost an explosion. So, for a calculation 
of spectra of evaporated particles with acceptable accuracy it is enough to know the value of a for an initial value of 
the PBH mass only. Taking this into account and having in mind the steepness of the PBH mass spectrum, we use 
the approximation 



a{m) = a{Af 



min \ 
BH ) 



(3.24) 



and just this value of a is meant in the expressions (3Tq)-(3.1f 
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TABLE I. Coefficients a,; and b, in Eq.(|3.25|) 



i 




u,d 


g 


s 


T 


c 


b 






t 


ai 


3.12 


9.36 


4.96 


9.36 


3.12 


9.36 


9.36 


1.86 


0.93 


9.36 


hi 


14.59 


14.13 


13.97 


13.89 


13.35 


13.32 


12.91 


11.87 


11.85 


11.39 



The very detailed calculation of the function a{m) was carried out in the works jllj,^. Here we use the simplified 
approach in which a{m) is represented by the dependence 

a = ao + ^ a, • 10^^ • 6(6. - logw{m)). (3.25) 

i 

Here, <d{x) is the step function. Coefficient ao gives the summary contribution of e+, e^, v and 7 and is equal to 
8.42 • 10^^ g'^sec-^ lllO- AU other coefficients are collected in the Table |. 

The coefficients hi determine the value of the PBH mass beginning from which particles of i-type can be evaporated. 
This value, M^^^, is obtained from the relations 

( GeV - kT'^H - ^ , (3.26) 



4 ,k(1)«6 



So, 



BH 



(3.27) 



(3.28) 



Writing Eqs. (|3.26 , ^.27 ) we took into account that the evaporation of massive particles becomes essential once the 
peak value in the energy distribution exceeds the particle rest mass. The peak value depends on the particle spin 
(see, e.g., 0): 



Epeak{s - 1/2) « AkT, 



BH 



'peak 



(,s ^ 1) 



6fcr, 



BH ■ 



(3.29) 



The values of masses of (u,d,s,c,b,g,T)-particles used in our calculations of are the same as in work [pT|, top quark 
mass was taken equal to 170GeV. The resulting function a{m) is shown on Fig.^ For comparison, the corresponding 
dependence from th e wor k (drawn using the Eq.7 of ||l^) is also shown. 

According to Eq.(3.14) the minimum value of PBH mass in the mass spectrum is about the horizon mass Mi. In 
turn, the value of Mi is connected with U, the moment of time just after reheating, from which the "normal" radiation 
era and the process of PBHs formation start, 



M, = ^ti 

8 tpi 



The corresponding reheating temperature is given, approximately, by the relation (see Eq.( 5.20 )) 

-1/2 



kTRH 



MeV . 



(3.30) 



(3.31) 



If kTnH values lie in the interval (10^ — 10^°) GeV (just for this interval we obtain the spectral index constraints in 
the present paper), the corresponding Mi values lie in the mass interval (10^^ — 10^^) g. So, as one can see, we really 
need for our aims the all information about a-parameter, that contains in Fig.^ including the highest values of a at 
m~ 10" g. 



The formula (3.18), as it stands, takes into account the contribution to the neutrino background solely from a 
direct process of the neutrino evaporation. If we suppose that particles evaporated by a black hole propagate freely, 
the calculation of other contributions to the neutrino background can be performed using our knowledge of particle 
physics lll|,0. 
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On Fig.|| the typical result of our calculation of instantaneous neutrino spectra from evaporating black hole is 
sho wn. The spectrum of the straightforward (direct) {ve + z^e) - emission is given by the Hawking function fniE, m), 
Eq.(|3.7D. The {v^ + t'e) spectrum arising from decays of (/x+ + ^^) evaporated directly is calculated by the formula 

f^''\E,m)^ J e{b^-\og,,m)9^fH{E^,m)^^:^^^^^^dE^ , (3.32) 

where c?n{^ /dE is the neutrino spectrum in a - decay, 5^ — 4. To evaluate the electron neutrino spectrum resulted 
from fragmentations of evaporated quarks, we used the simplest chain 

(m, d) — quarks — > tt — > ji — > Ve , (3.33) 

and the formula 

f'^'^>{E,m) = J e{bq-\og^om)gqfH[Eg,m) ■— ^^-^ ^ dEqdE^dE^, . (3.34) 

Here, gq = 24 for a summary contribution from u, u, d, d - quarks, dn'^/d^ is the g — s- tt fragmentation function, for 
which the simple parametrization was taken: 

and dnj^/ dE^ is the /i - spectrum in a decay tt — > /i + j/^. 

Analogous formulas are used for calculations of the neutrino spectrum from other channels of the neutrino produc- 
tion, for instanc e fro m decays of evaporated VF-bosons {W e-\- , W ^ ^ ^ i^e) ■ 

The formula ( 3.35| ) for the fragmentation function was suggested in the paper It leads to a simple multiplicity 



growth law of y E (su ch a law follows from a naive statistical model of an jet fragmentation ). It had been shown in 
ref. Hi that Eq.([3.35D fits rather weU the PETRA data and can be used instead of the more ingenious formula based 



on the leading logarithm approximation of QCD. In turn, it was shown recently |19|] using the programme HERWIG 
that the parametrizations of ref. being rather good at low_£, lead to an significant overestimation of the yield 



of large ^ final states. This fact as well as the simplicity of Eq.( 3.35 ) were decisive for our choice of the fragmentation 
function because, as we will show, we need just the upper estimate of the yield of quark fragmentations. 

In all calculations in the present paper we neglected the contribution to neutrino fluxes from neutron decays in the 
jets. These decays give the distinct bumps in instantaneous i^i>-spectra at i? ~ 1 MeV In general, neutrinos 

of such low energies have too small cross sections for neutrino-nucleus and neutrino-nucleon interactions (see sec.[v|) 
and are inessential for our main aim: the constraining of the spectral index (neutrinos with energies '--^ 10 — 100 MeV 
are mainly responsible for this). Besides, we will see at Sec.^ that at large values of one of the parameters of our 
approach (reheating temperature Tnu) these neutrinos are not noticeable in the summary neutrino background even 
aX E ^ 1 MeV because of strong redshift effects. 

In order to check the accuracy of our calculations of instantaneous neutrino spectra we calculated for one value of a 
black hole temperature the summary neutrino spectrum (excluding secondary Vt only). The result is shown on Fig.^ 
together with the corresponding curve from ref. [pl| . One can see from this figure that in the practically interesting 
region of neutrino energies {E > 100 MeV) our result differs from the Monte-Carlo result of ref. not more than 
on a factor of 2. 

The relative contribution of different channels to the total Ve spectrum strongly depends on the black hole tem- 
perature. One can see from Fig. 2, that at high temperature decays of massive particles evaporated by the black hole 
become rather important at high energy tail of the spectrum (if the corresponding branching ratios are not too small). 

The total instantaneous neutrino spectrum from a black hole evaporation is given by the sum 

f{E,m)=g,fH[E,m) + J^^'\E,m) + f(''\E,m) + ... , (3.36) 
and the total background neutrino spectrum is given by the same Eq.( 3.1^ ), except the change g^fniE, m) ,f{E, m). 



IV. NORMALIZATION OF THE PERTURBATION AMPLITUDE 

For numerical calculations of the neutrino background we need normalization of the dispersion ur entering the 
PBH mass spectrum formula. For this aim we must consider the time evolution of this dispersion. The connection of 
fii(i) with a power spectrum P{k,t) of the density perturbations is given by the expression 
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(4.1) 

oo 


Here, W{kR) is the smoothing window function. The second hne of this equation defines the (nonsmoothed) horizon 
crossing amphtude Snik, t). In general, it depends on the time and we need the initial amplitude Snik, U), whereas the 
normalization on COBE data (which had been performed, in particular, in the work |2l| ]) determines this amplitude 
at present time, t ^ Iq. Only in the particular case of a flat cosmological model with A the horizon crossing 
amplitude is time-independent and the normalization on COBE data gives Snik^ti) straightforwardly. 
Assuming a power law of primordial density perturbations, 

P{k,t)^k'', (4.2) 

and choosing for W{kR) the top-hat form, 

Wm = (4.3) 



one obtains from Eq.(4.1) 



4W = Cl^^SUkfut) EE ^^alikfut) . (4.4) 



In Eq.(4^), kfi is the comoving wave number, characterizing the perturbed region, 

kfi = \ (4.5) 

and C„ is a constant which depends on the spectral index n. For n ~ 1 one has C„ ~ 4.8 . The fluctuation power 
per logarithmic interval, or variance, at present time is given by the expression ||2^ 

^ ^ ,40) 

Here, As{k) is an amplitude of the scalar perturbations, k-dependence of which arises from a deviation from scale 
invariance (A| ~ fc""^). The factor g^iVlm) takes into account the suppression of the growth of density perturbations 
due to nonzero Oa (density perturbations s top growing once the universe becomes cosmological constant dominated) . 
The factor r2^„ in the denominator of Eq.(^^) arises because the potential fluctuations are reduced by ^m'- due to 
this, for fixed COBE normalization, the power spectrum of density perturbations must be, aside from other factors, 
enhanced by fi^^ 23 1. 

F or a rbitrary moment of time in matter era the expression for the power spectrum is evident generalization of 
Eq.(|!6|): 

|^P(M)^|A|,.)(»fe£^)\ (47) 

Here, the function g(a^ ^m) accounting for the growth of density perturbations is given by the expression p3| 



3 

5f2„i da f f da'\^ , 
5(a,l^™) = ^^/N^ da, (4.8) 



where r = Hot. In our case 

da ^ 1 



■ - + r^Afl^ (4.9) 
dr a 
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and 







(4.10) 



The normalization of this function is such that 



For the present time (a = Oo = 1) one has 



gia,l) = a . 



7(1, n,„) = g{flrn) , 



(4.11) 



(4.12) 



whe re g iflm) is the function entering Eq.( [4.6| ), as it must be. The horizon crossing amphtude Snik^t) defined by the 
Eq.(O) is determined by the relation 



[oHf 



Alik) 



g(a, ^m) 



Aiik) 



.g(^m) 



(4.13) 



One can see that the time dependence of 6H{k,t) arises just due to nonzero value of A. If A = 0, one has from 
Eq.(p^, using Eq.(Ell|), 



(4.14) 



Our normalization of perturbation amplitude on COBE data is based on two inputs. 

1. Connection between CMB fluctuations and scalar density fluctuations is described by the following relations 
(derived for a lowest-order reconstruction of the inflationary potential): 



5 = ^ =0.104/f (r!A)v42(fc,) , 
47r 



(4.15) 



4"^ (f^A) = 1.04 - 0.8217a + 2nl , 



(4.16) 



0.0 < r^A < 0.. 



(4.17) 



Here, C2 is the scalar contribution to the angular power spectrum of CMB temperature fluctuations for / = 2, 
fc* ~ ciqHq is the comoving wave number at the present horizon scale. The dependence fg*"^ on Qj^ is due to the 
integrated Sachs- Wolfe effect, i.e., due to the evolution of the potentials from last scattering surface till the present 
time. 

2. From 4-year COBE data one has (assuming that the scalar contribution dominates over the tensor on the 
large scales) the following results for small / values |2q]: 



2tt 



s\ 1/2 



(1.03 ±0.07) • 10" 



(4.18) 



n = 1.02 ±0.24 , 



From Eqs.(4.13-4.18) one obtains, finally, the normalization of the horizon crossing amplitude: 



fo(fc„to) = 4''°''''^-2-10-^ 



(4.19) 



(4.20) 



Asik^) =2-10 



-5 



(4.21) 
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All the dependence of 6H{k,t) on k is contained in an amplitude As{k) and, according to our assumption, Eq.(4.2), 
is very simple: 



As{k) ^ k' 



(4.22) 



Expressing the variable k through Alh, horizon mass at the moment when the scale A ~ crosses the Hubble radius, 
one obtains 



AsiMh) 



Mf^ " , radiation era , 



Mj^ " , matter era . 



(4.23) 



Writing Eqs.( 4.23 ) we took into account that the como ving mass density in radiation era decreases with time P7[ , 
while in matter era it is conserved. Now, using the Eq.(4.13) one can calculate the amplitude Snik) at the moment 
of matter - radiation equality: 



OH[Meq,teq) = — n As[Meq)- 



Here, M^q is the horizon mass at tf,q- Since a^q ^ 1, one has 



Hlq — HI—^ , g{aeq, flm) — ^eq ■ 



(4.24) 



(4.25) 



Therefore, according to Eq.(4.20), 

(5ff(Me„ie9) = AsiMh) = 2.0 • 10-^ 



1 ( Meg 



(4.26) 



Here, M^q is the present horizon mass. Using Eq.( [f.23| ) one obtains the time-independent amplitude SniMh) in 
radiation era: 



,5^(14) = 2- 10 



-5 



,/JI^)\MhoJ WqJ 



(4.27) 



The values of horizon masses M^q and Mho in Eq.(4.27) depend on Qm- 

Afeg = . 1.5 . io49^ ; Mho = n„,h-^ ■ 6 ■ lO^'^ g . (4.28) 
The expression for cj/f (Af;j) differs from Eq.( 4.27| ) only by the numerical coefficient (see Eq.(4.4)): 

<JH{Mh) = CnSniMh), (4.29) 



and Alh now is a function of k 



One can see from Eqs.(4.27 4.2q) that the horizon crossing amplitude aniMh) , normalized on COBE data, depends 



on f2„i only through the S^A-dependence of the fg''^ function: 



(JHiMh) 



I M- 



(4.30) 



The relations ( 4.27 ) and ( 4.30 ) are main results of this Section. The final expression for a nit) is, according to Eq.(4.4), 

'AW(t)\'/' 



Mfi{t) 



(4.31) 
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Here, Mfi{t) is the fluctuation mass, 



4 (a{t) 



3 



A^/;W=3^^-^) P{i) (4.32) 



and Mhor{t) is the horizon mass at a moment t. At the initial moment of time one has {MfioriU) = Mi , Mjiiti) = M) 

TV/A 2/3 



aniU) = aR{M) - ( ^ ) ^niMh) . (4.33) 



The expressions ( ^.llD of Sec-H can be obtained from Eqs.(|3o| , using the connection between M , Mi and Mh 



for radiation era (Eg. (2.1 



V. NEUTRINO BACKGROUND SPECTRA 

A precise calculation of the PBH neutrino background must include also a taking into account the neutrino absorp- 
tion during a travelling in the space. In the case of the photon background the most important absorption process is 
a pair production on neutral matter |l^ ] due to which the photons from PBHs evaporated earlier than z « 700 are 
absent today. In our case, the analog of an optical depth of the universe for the neutrino emitted at a redshift z and 
having today an energy E is given by the integra]|] 

z 

t{z,E)^ I a{E{l + z'))-n{z')^dz' . (5.1) 
J dz' 



Here, a{E) is the neutrino interaction cross section, n(z) is a number density of the target particles. 

Two processes are potentially " dangerous" : neutrino- nucleon inelastic scattering growing linearly with an energy, 
and annihilations with neutrinos of the relic background 

z/e + -/V — > + anything, (5-2) 
Ve + Ve{relic) ^^{f^ + ji). (5.3) 

i 

Here, fi are charged fermions (leptons and quarks). As is known, the relic neutrino background exists, with a Planck 
distribution, beginning from the epoch of neutrino decoupling {z k, 10^"). 

The cross section of the neutrino-nucleon inelastic scattering can be approximated (in the neutrino energy interval 
from E ~ IGeV to ~ lO^GeV) by the simple formula: 

a.AT (£;)« 0.5- 10-38 j^-^jcm^. (5.4) 

At i? > 10^ GeV the cross section grows with energy more slowly, so the use of Eq.( |5.4D in this region leads to the 
overestimation of the cross section value. But we will use this formula for the entire interval of neutrino energy, 
having in mind that practically we need only the upper estimate of this cross section (because, as we will see, the 
contribution of this channel to the total absorption factor is small). 

The number density of the target nucleons can be estimated using the existing constraints on fraction of critical 
density in baryons p8| |: 

0.004 < Os/i^ < 0.021. (5.5) 
From these inequalities and for h ~ 0.67 one obtains for the present nucleon number density : 



{l + zf 



= (4 • 10"® ^ 2 • 10"'^)cTO-3. (5.6) 



^ In formulas of Secs.^, ^ we use the convention h — c - 
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In numerical calculation we used the value 10 ^ for n^. The t^/at- function is given by the simple integral: 

z 

T^NiE, z) « 0.5 • 10-38 (^^^ . 10-7 . |(1 + zf^,dz' . (5.7) 



The expression ( ^.l[ ), as it stands, is valid only in a case when a target particle is at rest. In the case of an absorption 
through annihilations of evaporated neutrinos with relic neutrinos the current expression for the t^c, - function is (see, 
e.g., ref. ||) 

Z 

/I I 
((1 - cos9,^)a,o {E{1 + z),E,{l + z'),9,i,))n^{z') ■ —dz' . (5.8) 



The angular brackets indicate an average over the relic antineutrino energy distribution and over the angle O^o 
between the two colliding particles in the cosmic frame. 

The cross section of the vv - annihilation is expressed by the formula ]2^ , |30t 

a^,{E,E,,e,,) = ^ [N^f'jis) + N^ffis)] s, (5.9) 

s^2EEc{l^cos0^i^). (5.10) 

Here, s is a square of the summary energy of colliding particles in the center of mass system, the coefficients N^f / 
are the effec tive numbers of annihilation channels (for neutral and charged currents). 



Formula (5.9) is valid in the energy region before the Z boson pole, i.e., when 

s < M|. (5.11) 
As one can see from Eq.(p^), maximum values of s, for which we must know the cross section, are about E T~{l +z) 



2 



(where Tp is the present temperature of the relic antineutrino g as, z is the upper limit of the integral in Eq.(5.8), and 



E is the neutrino energy today). So, the cross section formula (5.9) can be used if 



M 



2 \ 1/2 / p \ -1/2 



^«(F%j -''"'{iv) ■ ("^' 

We will see from the numerical results that this condition is satisfied in all the practically important cases. 
The coefficients N^ff are given by the following expressions: 

^i// (s) = J2^^l- "^/) • (1 - sin^ Ow + Sqf sin" 0w) , (5.13) 

/ 



s 16 

N^j'j = e(- - "^e) • y sin^ Ow (5.14) 



Here, rif is the number of colours (1 for leptons, 3 for quarks ), t^f and qf are the third component of the weak 

isospin and the electric charge ( in units of the positron charge ), respectively, 9w is the electroweak mixing angle 
(sin^ 9w ~ 0.23), Q{x) is the step function. 
Using the relations 

{{l~cose,pfEp)^^pp — , (5.15) 
3 np 

{^~m})^ET9-m}, (5.16) 
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where po is the rehc antineutrino energy density, one obtains the final expression for t^c'- 

z 

T^.iE, z)- — - -pIE J{1 + z'f^,Nei} {ET^{1 + zT) dz' , (5.17) 



Neff^N^f'^ + N^/j. (5.18) 
Here, is the present electron antineutrino energy density, 

^ I . — (rP)4 . (5.19) 

Results of the calculations (with h = 0.67, fl^ — 0) of t^n and t^c are shown on Fig.^ for several values of neutrino 
energy. One can see, that the contribution of lyN channel to the total r is negligibly small everywhere. For typical 
neutrino energy ~ 100 MeV the absorption is essential beginning fro m z ^ 3 • 10^. A growth of the absorption factor 
with z is very fast, therefore the redshifts for which the condition ( ^.12| ) is violated are never really important in a 
calculation of the absorption. 

Our calculation of neutrino spectra from evaporating PBHs contains two parameters: a spectral index n and a time 
of an end of the inflation ti (which, by assumption, is a time when density fluctuations develop). We assume that at 
U the universe has as a result of the reheating the equilibrium temperature Trh- The connection of TfiH and ti is 
given by the standard model: 

+ n Qm -1/2 ^pi _ 0-24 ,^ 
t. = 0.301,. ^ ^ iT^H/lMeVr ' ^'-''^ 

(.9* ~ 100 is the number of the degrees of freedom in the early universe). 

On Fig.||the integrand of the background spectrum integral, considered in Sec. Ill is shown as a function of redshift z 



for several values of the parameter Tb,h- The inclusion of the absorption effects leads to an appearing of the additional 
factor in this integral: 

S{E)= [ e~^^'^-'^F{E,z)d\og^aiz + i) , (5.21) 



t{E,z) ^T^i,{E,z) + T^N{E,z) . (5.22) 

Each curve in Fig.^ has a strong cut off on some redshift value. This feature is connected with the existence in our 
model the minimum value of PBH mass. The PBH mass spectrum is steeply falling function of the mass, so the 
masses near minimum give a largest contribution to the neutrino background. The moment of their evaporation is, 
approximately, 

tev ~ Z ; (5.23) 

6a 

and the corresponding redshift is determined by the relation 

/t \'/' 

0e„ + 1 « {Zeq + 1) 7^ ■ (5-24) 

Larger masses evaporate at larger times and smaller redshifts. If, for instance, Tj^h = 10^° GeV, one has ti — 
0.24 • 10-26 s ; Af, = 7, 5 • lO^"^ g , ~ lO'^. So, at Trh = 10^° GeV the redshifts of order of 10^ give a largest 
contribution to the neutrino background (if there is no neutrino absorption; the absorption effects are essential just 
near z ~ 10^, and, as one can see from the figure, the actual z-dependences of the integrand have peak values at 
z ^ (2 3) • 10^, depending on a neutrino energy). 
The cut off value Zev strongly depends on Ti^h- 

^ ^ {MI"hT"^ - h'" - ■ (5.25) 
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Again, this estimate does not take into account the absorption effects. The real situation is seen on the figure. In 
particular, at Trh > 10^° GeV this cut off is the same for all Tjih values, as a result of a strong absorption at 
^ > 3 • 10^. 

On Figs.H,^ spectra of electron neutrino background arc shown separately for several channels of the neutrino 
production. One can see from Fig.^that 1) the contribution to the summary background from the direct Ve emission 
is, at TjiH }^ lO^GeV, absolutely dominant for E >1 MeV and 2) the relative importance of different channels changes 
with an increase of Trh- It is seen also that at Trh ^ lO^GeV the contribution of the quark fragmentation channel 
is very small at E ^ 1 MeV (< 5%), so the evident underestimation of this channel in our calculation, connected with 
the neglect of the contribution of heavy quark fragmentations, has (at large Tjih values) no particular importance. 

Situation with the contribution from quark fragmentations is slightly different at smaller Tj^h values (see FigJ^). 
Due to the flattening out of the curve for direct emission the relative contribution of the fragmentation channel 
at £^ ~ IMeV grows with an decrease of Trjj. The highest relative contribution is at Trh ~ 10^'^GeV. At this 
TnH the minimum value of mass in PBH mass spectrum (M]^jp) coincides with the threshold of quark evaporations 
{rrithr ~ 10^^-^g, according to Fig.|l|). At further decrease of Trh the relative contribution of quark fragmentations 
at ~ IMeV fastly decreases and, at some Trh value, it is strictly zero (when, due to correspondingly large M^jJ^ 
value, the age of the universe is not enough for "warming-up" the PBHs). 

It follows from Fig.0 that at neutrino energies larger than ^ (5-^- 10)MeV the contribution of fragmentation channels 
is negligibly small ( < 5%). As was mentioned above, the region of smaller neutrino energies is not important for the 
calculation of spectral index constraints, so some possible underestimation of the quark fragmentation contribution 
in this region (due to the simplified treatment of them in the present paper) is not sufficient. 

Some typical results of calculations of summary electron neutrino background spectra are shown on Figs.|,| The 
main features of such spectra have been revealed in previous works [ p^^lp2| : E~^ dependence above 100 MeV and 
the flattening out of the spectra at lower energies. The form of background spectra depends on the form of an initial 



PBH mass distribution. In refs. [Bl 32 1 the simple power mass spectrum of PBHs containing no cut off (minimum mass 



value) were used. The behavior of background spectra at low energies in this case is, aside from jet fragmentation 
contribution, ~ E~^ and turnover energy is ~ lOOMeV |^l[. In the case of PBH mass spectrum with the cut off 
the turnover energy depends on Trh [0 (it decreases as Trh increases, as is seen at F igja ). According to Figs.^,^ 
the turnover energy is about 1 MeV at Trh = 10^ GeV. Besides, it was shown in ref. p4[ that at low energies the 
photon background spectra are almost flat and, as is seen in Figs.H^, the same is also true in the case of the neutrino 
background spectra. 



VI. CONSTRAINTS ON THE SPECTRAL INDEX 



For an obtaining of the constraints on the spectral index we use three types of neutrino experiments. 

1. Radiochemical experiments for the detection of solar neutrinos. There are data on solar neutrino fluxes from the 
famous Davis experiment | |3^ and the Ga — Ge experiment ]3^ . The basic process used in radiochemical experiments 
is a charged-current neutrino capture reaction 



v, + A''z^e-+A%-i. (6.1) 

In general, the cross section for the neutrino absorption via a bound-bound transition can be calculated using the 
approximate formula 



cr = < 




PeEe , E < lOOMeV 

(6.2) 



. const. , E > lOOMeV. 

Here, pe and E^ are the momentum and the energy of an electron produced in a neutrino capture reaction. Neutrino 
and electron energies are connected by the relation 

E, = E- Ethr + me . (6.3) 

Stopping of the rise of a neutrino capture cross section (for a bound-bound transition) with neutrino energy is due to 
influence of a weak nuclear formfactor of the transition. The approximate value of the turnover energy in Eq.(3.2) is 
determined from the relations 

E ~ Qmax ~ 7^ ~ TOtj- ■ (6.4) 
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Here, qmax is a maximum momentum transfer from a neutrino to the nucleus (in a model with step-like nuclear 
formfactor), Ra is a size of the nucleus. 

In the case of the Davis experiment the neutrino capture reaction is 

CI Ar + e- . (6.5) 

Here we took into account the super- allowed transition only, i.e., the transition to the isotopic analog state in ^^Ar, 
for which 

(1)2 ^ 3 ^ (cr)2 ^ 0.2 , Ethr ~ 5MeV. (6.6) 

This transition gives dominant contribution to the effect in the detector, in spite of the rather high energy threshold, 
due to the large value of the Fermi strength (l)^. 
In the case of the Ga — Ge reaction, 

f Ga Ge + e- , (6.7) 

the main contribution gives the ground state - ground state transition {Ethr — 0.242MeV). In ^^Ge also there is an 
isotopic analog state, but its energy is too high, ~ 9 MeV (it is slightly higher than the particle emission thr esho ld, 
so this transition is irrelevant). The cross section of the ground state - ground state transition in the reaction ( |6.7| ) is 
given by the expression (for E < 100 MeV) 

^ 2Ja + 1 6163.4 

Here, Ja and Jb are an gula r momenta of ^^Ga and ^^Ge, respectively, ft value is determined by the half- life time of 
^^Ge. According to ref. |^ one has 

Ja^I , Jb = 1 , lg/t = 4.37 . (6.9) 

The number of target atoms N'^ in the detectors is about 2.2 • 10'^° for the Gl — Ar experiment |Q and 10^^ 
for the Ga — Ge experiments |^^. The average statistics is '--^ 1.5 captures/day {Gl — Ar) and 1 capture/day 
(Ga — Ge). The constraints are calculated using the relation 

Ait ■ ■IQr' ■ j S{E)a{E)dE <\. (6.10) 

2. The experiment on a search of an antineutrino flux from the Sun In some theoretical schemes (e.g., in 

the model of a spin - flavor precession in a magnetic field) the Sun can emit rather large flux of antineutrinos. LSD 



experiment |36| sets the upper limit on this flux, <f>£,/$i^ < 1.7%. In this experiment the neutrino detection is carried 



out using the reaction 

i>e+p^n + e~^ . (6.11) 
The number of target protons is ~ 8.6 • 10^* per 1 ton of the scintillation detector, and the obtained upper limit is 



0.28 antineutrino events per year per ton The cross section of the reaction (6.11) is well known (see, e.g., [p7[). 
In particular, in the low energy region {E ^ mp) it is given by the formula 



a{^ + p + e+ + n) = ^\l + 3[^] ] p,E„ (6.12) 




Ee = E - {nin - nip). (6.13) 
The threshold of this reaction is given by the relation 

Ethr = m.n -mp + nie^ 1.8 MeV. (6.14) 



The cross section of the reaction ( ^.11 ) grows with the neutrino energy up to i? ~ 2 GeV, and is a constant 



(~ 0.5 • 10 cm^) at larger energies. The product of this cross section and a PBH antineutrino background spectrum 
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has the broad maximum at E lOOMeV. The constraint is determined from the condition that the calculated effect 
in the LSD detector is smaller than the upper limit obtained in ref. p6|] . 

3. The Kamiokande experiment on a detection of atmospheric electron neutrinos fs^ . In this experiment the 
electrons arising in the reaction 

i/f™ + n^39 + e- (6.15) 

in the large water Cherenkov detector were detected and, moreover, their energy spectrum was measured. This 
spectrum has a maximum at the energy about SOOMeT^. The spectrum of the atmospheric electron neutrinos is 
calculated(see, e.g., [^) with a very large accuracy (assuming an absence of the neutrino oscillations) and the 
experimentally measured electron spectrum coincides, more or less, with the theoretical prediction. The observed 
electron excess at ~ IQQMeV (which is a possible consequence of the oscillations) is not too large. We use the 
following condition for an obtaining the our constraint : the absolute differential intensity of the PBH neutrino 
background at the neutrino energy E ~ Q.iGeV cannot exceed the theoretical differential intensity of the atmospheric 
electron neutrinos at the same energy (otherwise the total number of electrons in the detector and the energy spectrum 
of these electrons must be different from the observed ones) . 



VII. RESULTS AND DISCUSSIONS 



Fig .|lO|,|ll| shows our results for the spectral index constraints. It is seen that the best constraints are obtained using 
the Kamiokande atmospheric neutrino data and the LSD upper limit on an antineutrino flux from the Sun. 

It appears that the spectral index constraints obtained in the present work are not stronger than the corresponding 
constraint following from COBE data (see Eq.( 4.19| )). One should remember, however, that our approach directly 
constrains the perturbation amplitude aniMh) rather than the spectral index. Comoving length scales of pertur- 
bations, for which we constrain the amplitude, are very small in comparison with the size of present horizon. For 
example, if Trh ~ 3 • 10^ GeV the minimum PBH mass value is about 10^^ g. In this case the time of PBH formation 
is 



U 



pi 



pi 



10' 



(7.1) 



the corresponding scale factor is, according to Eq.(3.21), about^^ 10 , and for the comoving length scale one obtains 



a 



10"^ pc. 



(7.2) 



It is clear that the normalization on COBE data has sense only if it is assumed that the spectral index is independent 
on scale. The curves of Fig.|l^,|ri| were obtained using just this assumption; it is very easy to calculate from them the 
corresponding constraints on auiMh). 

For a demonstration of the sensitivity of these results to a chosen f^A-value we calculate the constraints for two 
extreme J^A-values, for which the parametrization for (Eq.( 4.1^ )) is still valid. It is clear that some weakening of 
the constraint with an increase of JIa is due to the corresponding decrease of values of the function <t//(M/j) entering 
the initial PBH mass spectrum (see Eq.( 2.20| )). This decrease of aH{Mh) takes place solely due to the dependence of 



4°^ on f^A (see Eq.(|l6|)). 



The behavior of the constraint curve n{TjiH) is sharply different from that was obtained in ref. |15[ | (authors of ref. 
JlSl used the initial PBH mass function following from the near critical collapse scenario with a domination of the 
earliest epoch of PBHs formation, and diffuse extragalactic photon data). 

The slight bend at T^h ^ 10^"GeV of the constraint curves (Fig.|l0|) is an effect of a neutrino absorption in the 
space. Constraint based on the atmospheric neutrino experiment is especially sensitive to the absorption because 
neutrinos of larger energies (in average) are responsible for it. 

Flattening out of the con straint curve n iTj^n) with increase of Tj^h is connected with the dependence M'^^{Tjih)- 
According to the relations ( 3.14 , 3^30| , 5.20 ), the minimum value of PBH mass in the initial mass spectrum is inversely 
proportional to T^j^. More exactly, the relation between M]^Jp and Thh is 



M^if(g) 



7.2 • 10 



30 



(TRH/lGeV)^ 



(7.3) 
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At small TijH the minimum value of PBH mass in the initial PBH mass spectrum is large and the neutrino background 
is dominated by the evaporations at recent epochs. In opposite, at large Trh, when Mg^ is small, the neutrino 
background is dominated by the evaporations at earlier times. Thus, using the relation 

~ 2.4 • l0^n,r^h^ (7.4) 

and the curves, corresponding to T^h = lO^GeV in FigJ|, one can see that in cases when Tj^h > lO^GcV the neutrino 
background is dominated by PBH evaporations during the radiation era. 
Introducing, as usual [Q, the function M*(M]^|p) by the formula 

M* {M^^^) = [MMTH)tA , (7.5) 

where is the age of the universe, one can easily see that at low Tjih (before the flattening out of our n(TiiH) curve) 
the inequahty M]^ip > M*{M^'i}') takes place while at large Trh one has M^}p < Af*(MfJp). We remind that M* 
is the mass of a black hole, formed in the early universe and evaporating today. The turning-point (i.e., the point of a 
change of the slope of the n{TjiH) curve) corresponds to the equality M]^Jp — M*{M'^^) and in the case SIa = it is 
at TfiH ~ lO^GeV. The shift of this turning-point with an increase of 0,^ is explained by the corresponding increase 
of an age of the universe: 



to = = In ^ " . (7.6) 

One can compare our spectral index constraints with the corresponding results of ref. |p^ , where the same initial 
mass spectrum of PBHs had been used. Some difference in resulting constraints (An ^ 0.02) may be connected simply 
with the fact that authors of ref. use slightly different formula for aH{Mh), namely, 

aif(M,.) = (Tff(A4,o) fi^) ' . (7.7) 



Mh . 

In other respects the constraints are quite similar although in ref. |l4[| they were based on diffuse extragalactic photon 
background data. 

One should note that usually the calculations of these constraints are accompanied by the calculation of the bounds 
based on requirement that the energy density in PBHs does not overdose the universe at any epoch {^bh < !)■ For 
a setting of such bounds one must consider the cosmological evolution of the system PBHs -I- radiation. We intend 
to carry out these calculations in a separate paper. Estimates show that at large Tfin values {Thh ^ 10^'^GeV) the 
constraints based on Qbh < 1 are stronger than those based on the neutrino experiments (it is the reason why we 
did not explore the region Tj^h > 10^'^GcV in the present paper). 

The neutrino background spectra obtained in the present paper may be useful somewhere else, irrespective of the 
spectral index constraint problem. The main feature of such spectra a.t E > IMeV and at Tpjf > 10^ — 10^" GeV is 
the very small contribution (< 5%) of all secondary channels of the neutrino production at evaporation (except the 
decay of directly evaporated muons) (see Fig.^. There are two reasons for this: the steepness of the initial PBH mass 
spectrum and, most essentially, the existence in our model the cut off in this spectrum, i.e., the minimum PBH mass 
value. Evaporation of PBHs having masses near M"^^ gives the dominant contribution in the background, so there 
is a strong maximum in z distributions at low energies (Figj^). For example, at Tfjji = 10^'^GeV the peak value of 
2; is ~ 10^. Neutrino energy ~ IMeV corresponds in this case to energy at evaporation ~ ITeV. For an evaporation 
of such high energy neutrinos the PBH mass must be ~ 10^°g in the case of direct neutrinos and much smaller in 
the case of secondary neutrinos (because their energies are much lower in average (see FigJ|)). At the same time, 
the minimum value of PBH mass in the mass spectrum at Trh = 10"'^''GeV is ~ lO^^g. So, it should be clear that, 
in our example, the secondar y neu trinos with E ~ IMeV can appear only after the essential evolution of the initial 
PBH mass spectrum (see Eq.( ^.17 ) and the detailed discussion of this evolution in ref. |^^). This evolution leads to 



an increase of the total mass interval and to the corresponding decrease of a PBH number density (per unit mass 
interval) and, as a result, to a relative smallness of the secondary neutrino flux. 

One can compare our results for the neutrino background spectra with those of ref. | |3l| , |3^ . According to the 
conclusions of ref. |^^, the secondary neutrino fluxes are larger relatively and even dominant a,i E < lOMeV. Similar 
results were obtained in earlier work In our model the dominance of secondary neutrino fluxes at < 10 MeV 
takes place at some values of Trh only (see Fig.^ . 

The results for neutrino diffuse background spectra and the constraints on the spectral index obtained here do 
not take into account the possible existence of the photospheres [p2|-p4| around a black hole. In studies of photon 
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background spectra two kinds of photospheres should be considered: QCD photosphere (quark-gluon plasma arising 
from interactions of quarks and gluons before hadronization) and QED photosphere (e"*" — e~ — 7 plasma, in which 
directly evaporated photons are processed). As for the latter, in our case one should consider instead the "electroweak 
(or neutrino) photosphere" |42| which may form only at extremely high black hole temperatures and, probably, has 
no influence on the background spectra. QCD photosphere may be more important for the problem. Thus, the 
recent calculations |Q] of instantaneous photon spectra from the model QCD photosphere show that some increase 
of photon fluxes (in comparison with the usual case when QCD photosphere does not form and photons appear simply 
through the direct quark fragmentations with subsequent Tr'^-decays) definitely takes place. For example, according to 
these calculations, at Tbh = 50GeV the peak value of instantaneous photon spectrum increases on a factor ^ 4. If it 
is true and if this increase is the same at all black hole temperatures, the contribution to neutrino background spectra 
from the quark fragmentations obtained in our calculations must be increased approximately on the same factor. 
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FIG. 1. The function a{m) counting the degrees of freedom of the PBH radiation (m is an instantaneous value of the PBH 
mass). SoUd line is the present result, dashed line is drawn using Eq.(7) of Ref.[17]. 
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FIG. 2. Instantaneous electron neutrino spectra from the evaporation of PBH with Tbh = ITeV, for several channels of the 
neutrino production. 
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FIG. 3. Instantaneous neutrino spectra from a black hole with a temperature of 100 GeV. 



23 



N 

w 

e 

o 

o 



4 - 



-2 - 



-4 



-6 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



^A=0 




lO^GeV 



10 ^GeV 



10 ^GeV 



VV 



10"GeV 



10 ^GeV 



10 ^GeV 



VN 



I .... I .... I .... I 



J I l_l I l_l I I L_I I l_l I 1_1_ 



4 4.5 5 5.5 6 6.5 7 7.5 8 8.5 9 



log.,(z+l) 



'10 

FIG. 4. Absorption factor for two channels of neutrino absorption and for several values of neutrino energy, as a function of 
0, the redshift corresponding to a moment of the neutrino emission. 



24 



-4 
-6 
-8 
-10 
-12 
-14 
-16 
-18 
-20 
-22 
-24 





* 










10"^ GeV 












y' 10"^ GeV 
















L Tjyj=10^-^GeV n=1.25 

1 .... 1 . 





> 
o 

C5 



N 

w 



01- 
O 



2 

logi„(z+l) 




Tjy,=10 GeV n=1.25 
■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ 



4 5 
logi„(z+l) 



6 
4 
2 

-2 
-4 
-6 
-8 
-10 
-12 
-14 
-16 
-18 
-20 



1 1 1 1 1 


■ ■ 1 ■ ■ 1 1 ■ 






^ " 


\ 








- 10-^ 


GeV y \ 


















10"^ GeV 




* 

* > 






* ^ 












~ * y 

- ' / 

\ *y 






- T^=10'°GeV 
. . . 1 


n=1.25 

. . 1 ... 1 . 





> 



u 



w 



01 

o 



6 8 
log,.(z+l) 




T^=10 GeV n=1.25 
J , , I , , , I .. 



6 8 
log^(,(z+l) 



FIG. 5. Redshift dependence of the integrand of the expression (5.21) for a neutrino background spectrum (with JIa 
for two values of the neutrino energy and for several values of the parameter Trh ■ 



25 




log,. (E/lGeV) 



_ 8 

ca 
m 

H 

N 

1' -2 



-4 



Tj^=10 "GeV n=1.25 




(b) 



direct V 



direct l^-^v^ 
W^^Vg and W*^n,^Vg 

. I .... I .... I .... I ... . 
-3 -2 -1 



log,. (E/lGeV) 
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FIG. 8. Electron neutrino background spectra from PBHs, calculated for several values of reheating temperature. Dashed 
curve shows the spectrum for TjiH = 10^°GeV, calculated without taking into account the neutrino absorption. 
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FIG. 9. Electron neutrino background spectra from PBHs, calculated for several values of the spectral index. Dashed curve 
shows the theoretical atmospheric neutrino spectrum at Kamiokande site [39] (averaged over all directions). 
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FIG. 10. Constraints on the spectral index n as a function of the reheating temperature Trh, from three types of the neutrino 
experiments. 



30 



gl.29 



1.28 - 



1.27 



1.26 



1.25 - 



1.24 - 



1.23 



-1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



a^=o . 8 



_i I I I I I I I I I I I I I I I I I I I I I I i_ 



7.5 



8 



8.5 



9.5 



10 



login (T.H/lGeV) 



'10 ^ RH' 

FIG. 11. Constraints on the spectral index n as a function of the reheating temperature Trh, based on the solar i>e experiment, 
for two values of Qa- 



31 



